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Extensive LAPW frozen phonon calculations were performed in order to 

On 

understand the origin of ferroelectricity in LiTa03 and L1NDO3. Displacement 
of the Li atoms alone results in an anharmonic single well, whereas displace- 
ments of oxygen and lithium together result in deep double wells, much deeper 
than the transition temperatures, Tq- This is contrary to current theories 

x 

5_i ■ which model the underlying potential as a triple well potential for the lithium 

atoms. Our results support an order-disorder model for the oxygen atoms as 
the driving mechanism for the ferroelectric instability. Oxygen displacements 
alone against the transition metal atoms result in shallower double wells as a 
result of oxygen-lithium overlap so that the lithium and oxygen displacements 
are strongly coupled . We find large hybridization between the oxygens and 
the transition metal atoms. Thus ferroelectricity in the Li(Nb,Ta)03 system 
is similar in origin to ferroelectricity in the perovskites. We also find that 
the electronic structures of LiTa03 and LiNb03 are very similar and hardly 
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change during the phase transition. 
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I. INTRODUCTION 



The origin of ferroelectricity in the two well known ferroelectric systems, LiNbOs and 
LiTaC>3, has been subject to intense study since the discovery of LiNbOs in 19490. They have 
many applications in optical, electro-optical and piezoelectric devices, but the fundamental 
physics that leads to their ferroelectric behavior have not been studied. Their transition 
temperatures, which are among the highest known ferroelectric transition temperatures, are 
quite different, 1480 K for LiNbOs and 950 K for LiTaOs. The electronic origin of their 
different Tq is a mystery since Nb 5+ and Ta 5+ behave very similarly, and structurely these 
materials are almost identical. The origin of their ferroelectric instability as well as their 
different transition temperatures is investigated. 

Both materials undergo only one structural phase transition. The paraelectric structure 
has a 10 atom unit cell and the average structure belongs to the R3c space group, The atomic 
arrangement consists of oxygen octahedra sharing faces along the polar trigonal axis. The 
transition metal atoms occupy the centers of oxygen octahedra, and the average Li atom 
position lies on the face between two adjacent oxygen octahedra (Fig. |]a). The ferroelectric 
structure is rhombohedral, and belongs to the space group R3c. The transition metal atom 
is displaced from the center of the oxygen octahedra along the trigonal axis. The next 
oxygen octahedron along this axis is empty and the adjacent octahedron has a Li atom 
ferroelectrically displaced from the oxygen face in the spontaneous polarization P s direction 
(Fig. |l|b). Glass in 1968 and later Johnston et al. have determined, using dielectric and 
thermal measurements, that the phase transformation in these systems is continuous!'!. 

Whether the transition is displacive or order-disorder has been much discussed and con- 
fusion abounds. A displacive phase transition is one where the local potential in the mean 
field of the rest of the crystal has a single minima, and is characterized by a temperature 
dependent optic mode approaching zero as the temperature reaches Tc- Temperature de- 
pendence measurements of Raman, Rayleigh scattering)! and infrared reflectivitylH! show 
soft mode behavior for one polar Ai (TO) optic mode in the ferroelectric phase. This soft 
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mode crosses many E modes (whose eigenvectors give ionic displacements perpendicular to 
the polar axis) , and thus this specific mode is difficult to trace in detail. Tomeno and 
Matsumura@ measured the dielectric constants of LiTaC>3 and found a large Curie constant, 
and interpreted their results as indicative of a displacive transition. 

Evidence for the transition having an order-disorder character came from Penna et aZ.§H0 5 
who observed no mode softening for LiTa03 for the Ai (TO) mode, followed by Chowdhury 
et a/.0 who performed neutron scattering experiments on LiNbOa, and also failed to observe 
any softening of the Ai mode. Okamoto et a/.0 used Raman scattering to study LiNbOs 
between room temperature and 1225K and saw two of the three Ai modes (which are, at 
room temperature, TO, LO and TO at 250, 270 and 274 cm -1 , respectively), and observed 
anomalous behavior for one of them (at 274 cm -1 ) as the temperature reached T c . They 
noted that the decrease in the peak frequency was mostly due to the rapid increase in 
damping as the linewidth had a divergent form in temperature, whereas the quasiharmonic 
frequency remained almost temperature- independent. Their conclusion was therefore that 
LiNbOa does not exhibit a typical displacive transition, but rather resembles an order- 
disorder system. Zhang et a/.0 reported these measurements on LiTaOs, and found the same 
kind of behavior. In an order- disorder phase transition the local potential is characterized by 
a double (or more) well, with the thermal energy kTc much smaller then the well depth and 
no soft phonon mode exist since phonons now oscillate within each well and the wells remain 
essentially unchanged throughout the phase transition. These transitions are characterized 
by a diffusive soft mode which is not a phonon but represent large amplitude thermal hopping 
between the wells. At T above Tq the crystal is non-polar in a thermally averaged sense. 

Jayaraman et a/.lll argued for an order-disorder type transition because they saw little 
pressure dependence of the Raman mode; they also emphasize the difference from ferroelec- 
trie perovskites, which show a strong pressure dependence. RaptisB measured aud analyzed 
Raman modes of LiTaOs between room temperature and 1200K and observed softening of 
this Ai mode (along with others) to a certain degree. However, the decrease was charac- 
terized with an order-disorder model. Catchen et a/.0 used perturbed-angular-correlation 



(PAC) spectroscopy to measure nuclei-electric-quadrupole interactions at the Li sites over a 
temperature range of 295-1 100K, and Cheng et al. fE! studied inelastic neutron scattering 
from room temperature up to 800°C (1100 K) in LiTaOs; both failed to observe mode soft- 
ening, therefore not supporting the displacive picture for the phase transition. Tezuka et 
aZ.El used hyper-Raman and Raman spectra of LiTa03 between 14 and 1200 K. No evidence 
was found for the softening of the A\ mode; however, a strong Debye-type relaxational mode 
was found in the two phases, suggesting an order- disorder type transition. They interpreted 
the anomalous line shape of an A\ mode in terms of coupling with relaxational modes. 

Most ferroelectric systems are thought of as exhibiting displacive behavior far from 
the transition temperature region and order- disorder characteristics near Tq. In the 
Li(Nb,Ta)0 3 systems this conclusion is supported by a number of studies&M 

In other ferroelectric oxides, like the perovskites KTaC>3 and KNbOs, the mechanism 
behind the phase transition has also been debated. Evidence for the transition being of the 
displacive type are presented by Nunes et a/.0 and Samara@, while Comes et a/.@ and PAC 
experiments by Dougherty et a/.il point to it being of the order-disorder type. Sokoloff et 
have studied Raman scattering of KNbC>3 and BaTiC>3 and discovered central peaks 
which have line shapes and thermal dependence characteristics of Debye relaxation modes 
as well as symmetry properties consistent with the eight-site model. A theoretical study by 
Edwardson of KNbOsil, using interacting polarizable ions in static and dynamic simulations 
found a mixture of order-disorder and displacive type behavior. Postnikov et a/.@0 carried 
out an LMTO study examining the total energy of KNbC>3 in the tetragonal and rhom- 
bohedral phases. These calculations found that displacing the Nb atom along the <100> 
direction from its tetragonal cell position corresponds to a saddle point on the total energy 
surface. This becomes a minima when the lattice strain is included. They have also carried 
out total energy calculations for KTaOs which found no ferroelectric instability for the cal- 
culated volume. They were able to induce a phase transition by applying negative pressure 
(expanding the lattice). An XAFS study of KTN by Hanske-Petitpierre et a/.0 found that 
the ferroelectric transition is not displacive and involves orientational order- disorder transi- 



tion of the Nb atom. A recent first principles investigation of eight perovskitesil'il suggests 
that in materials like KNDO3 and BaTi03, which are rhombohedral at T=0, the sequence 
of successive transitions is explained via the eight site model, where the order parameter in 
the paraelectric phase fluctuates between the eight minima in the [111] directions. These 
sites are minima at the cubic phase, before the development of the strain. 

The theories developed thus far for the LiTaC>3 and LiNbC>3 systems usually are based on 
the Lines modefi M. E. Lines applied his effective-field theory to LiTaOs, and parametrized 
it as a displacive ferroelectric due to the data available at the time by Johnston and 
Kaminowi, and assumed a triple well potential for the Li atoms. Abrahams et a/S per- 
formed neutron scattering of LiTaC>3 between room temperature and 940 K, and discovered 
that, above Tq, the lithium atom positions in LiTaC>3 become disordered and hop among 
the centrosymmetric position and sites at ± 0.37 A along the optic axis. Similar measure- 
ments for LiNbC>3 show the same kind of behaviorll. The neutron scattering data is the 
cornerstone behind all theories modeling this ferroelectric transition as an order- disorder 
mechanism with the Li ions hopping among the centrosymmetric sites and the adjacent oc- 
tahedral sites. This approach was adopted by Birnie&il, who modeled the Li hopping as a 
Frenkel defect, and later by Bakker et a/.0, who used this data in addition to the triple well 
Lines model as a basis for a quantum-mechanical description of the phase transformation in 
LiTaC>3. Bakker et al. predicted and observed^ a 32 cm -1 excitation which they ascribed 
to Li motions between the central and lowest wells. This excitation frequency has not been 
observed in other studies, howeveiil. 



II. METHOD 



The Kohn-Sham equations^S are solved self-consistently using the full potential Lin- 
earized Augmented Plane Wave (LAPW) method^, where the electronic many-body 
exchange- correlation interactions are described by the local density approximation (LDA) 
using Hedin-Lundqvist parametrization^. There are no shape approximations for the charge 
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density or the potential. This method has proved predictive in many previous studies. Ex- 
amples include the prediction of a high pressure phase transformation in silicaEzI, studies 
of iron at high pressurest3l±3 studies of AI2O 30, MgO and CaO0, MgSi0 3 i and high tem- 
perature superconductors^. This method was previously applied successfully to perovskite 
ferroelectrics like BaTi0 3 , PbTi0 3 @ and KNbojl. 

We use the LAPW+LO methodH which uses a mixed basis consisting of the LAPW 
basis plus extra localized orbitals inside the muffin tin spheres. The extra local orbitals 
remove a Li s ghost state and relax the valence states. It also allows the use of a single 
energy window. Local orbitals included s for Li, s and p for O and s, p and d for the Nb 
and Ta atoms. Other details of the calculations include a muffin tin size of 1.6 Bohr for the 
Li and 2.0 Bohr for the Nb and Ta atoms. The oxygen's muffin tin radius was 1.6 Bohr for 
oxygen-Nb separation up to 1.882A (3.556 Bohr). The corresponding distance in LiTa0 3 is 
O-Ta up to 1.8845A (3.562 Bohr). At this point the oxygen muffin tin radius was decreased 
to 1.552 Bohr. For smaller separations the muffin tin radius was 1.506 Bohr. In order to 
be able to compare the energies calculated using different muffin tin radii, we have repeated 
calculations with the three sets of muffin tin radii to find the energy shift due to this change 
in the muffin tin radii, and have assumed that this small shift (~ 7 mRy) is constant for 
small displacements of atoms. 

A 4x4x4 special k-point mesh was used which generates a total of 10 k-points in the 
irreducible zone. To test energy convergence, the energies at the symmetric and experimen- 
tal structures of both LiTa0 3 and LiNb0 3 were also calculated with a 6x6x6 mesh which 
generates 28 k-points in the irreducible zone, and these energies are shown in Table |. The 
change in energy difference for the two k-point sets between the experimental and the sym- 
metric configurations is 0.069 mRy for LiNb0 3 and 0.3 mRy for LiTa0 3 , demonstrating 
convergence. 

The RK max parameter was set to 7.0 which gives approximately 1150 basis functions for 
the LiTa0 3 calculations and 1050 functions for LiNb0 3 . The core states were calculated 
fully relativistically and the valence states semi-relativistically. For each Ta atom, the states 
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up to 4f were included in the core, and as a result 0.588 electrons extended beyond the 
muffin tin sphere. For each Nb atom states up to 4s were included in the core and only 
0.07 core electrons extended beyond the sphere; core electrons that spill out of the muffin 
tins see an extrapolated spherical core potential. Also, in LiNb03 the Li atom s states were 
included as bands, whereas in LiTa0 3 they were treated as core states. 

III. RESULTS AND DISCUSSION 
A. Energetics 

We have calculated the potential energy surfaces along the experimental soft mode coor- 
dinate. To test the sensitivity to the different lattice parameters, the total energy of LiNb0 3 
was calculated in the ferroelectric configuration using both the LiNb03 lattice parame- 
ters (a^=5.14829A and c H =13.863lA )H, and the LiTa0 3 lattice parameters (a^=5.15428 
A and c#=13. 78351 A )0, a difference of 0.75% in the c/a ratio. The effect of this strain on 
the total energies was almost negligible; slightly less then 1 mRy or 5.5% of the well depth. 
This is in contrast to the case of the perovskites; e.g. PbTi0 3 , where a strong dependence of 
the total energy on the tetragonal strain was observed, and the energy decreases markedly, 
about 35% of the well depth for the experimental 6% c/a stram@. Also the total energy of 
LiNb0 3 using both the experimental LiNb0 3 atomic positions (Li at (0.2829,0.2829,0.2829), 
Nb at (0,0,0) and Oxygen at (0.1139,0.3601,-0.2799) )0 and the LiTa0 3 atomic positions (Li 
at (0.279,0.279,0.279), Ta at (0,0,0) and Oxygen at (0.1188,0.3622,-0.2749)1) were calcu- 
lated. The resulting wells are less then 1 mRy different, the LiNb0 3 atomic positions yielding 
the deeper well. For the purpose of comparison, except for a few more points which yielded 
the same results (energy differences of less then 1 mRy), all points were calculated using 
the experimental lattice parameters and positions of LiTa0 3 to facilitate comparison of the 
effects of chemistry on ferroelectric behavior and electronic structure. Table [I] summarizes 
the results for LiTa0 3 and LiNb0 3 . The first column refers to the displaced atoms, and to 
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the amount of displacement as a fraction of the paraelectric to experimental ferroelectric 
normal mode amplitude. Note that one distortion of LiTaC>3 published earlierll was not 
along the soft mode coordinate; the present results correct this error. 

Figure ^|a shows the potential energy surfaces of LiTaC>3 with respect to displacements 
of Li only (upper curves), oxygen only (middle, shallow double wells) and Li+O (lower 
curves). Figure |2|b shows the same picture for LiNbC>3. The lithium displacements along 
the soft mode coordinate result in a single anharmonic well with low curvature. Displacing 
only the oxygens against the transition metal atoms results in shallow double wells, and 
the deep double wells are the result of the Li+O displacements along the experimental 
ferroelectric coordinate. 

The wells resulting from the oxygen and lithium displacements have well depths of 17.3 
mRy (2739 K) and 18.3 mRy (2858 K) for LiTa0 3 and LiNb0 3 , respectively. Both wells 
are much deeper than the experimental transition temperatures, which is consistent with an 
order-disorder character for the phase transition. 

The energy was fit to a fourth order polynomial in normal mode amplitude, Q = 
yJ2i rriiiif. The Schrodinger equation was solved numerically to obtain the eigenstates 
assuming one- dimensional non-interacting anharmonic oscillators along the soft-mode co- 
ordinate. Figure [| shows the energy surface as a function of the normal mode amplitude, Q, 
and the energy levels. We can see that the two wells have a different shape due to the factor 
of about two in mass of Nb and Ta, and thus the Nb approximately displaces twice as much 
as the Ta relative to the center of mass, which results in a different normal mode amplitude 
Q. The energy difference between the ground and lowest excited state gives a frequency 
for LiTaC>3 of 270 cm -1 in fairly good agreement with the experimental Raman frequency 
of 201-225 cm -1 considering the one dimensional noninteracting oscillators approximation. 
For LiNbC>3, the calculated frequency is 250 cm -1 , in excellent agreement with experimental 
data of about 250-275 cm -1 . 

These results indicate that these structural phase transitions are not dominated energet- 
ically by the displacements of the lithium alone. The potential energy surfaces show that 
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the deep double wells are the result of the coupled motion of lithiums and oxygens. Dis- 
placement of the lithiums alone hardly changes the energy of the system. This is in contrast 
to current theories which model the displacement of the lithiums as the driving mechanism 
for the ferroelectric instability. 

In order to understand the oxygen-lithium coupling, we calculated the dynamical matri- 
ces for the LAPW and Madelung energies (assuming fully charged ions) for both materials. 
The LAPW and Madelung energies were fitted to a 4th order polynomial surfaces in the 
normal mode coordinates of the lithium and oxygen (Qli and Qo)- The second derivatives 
of these energy surfaces at zero displacements are the coefficients of the dynamical matrices. 



Table [TTJ shows the coefficients of the fit for the total energy of LiNbC>3 and LiTaC^; all 
the coefficients for LiNb03 and LiTaC>3 are well constrained except the coefficient of Q 2 Li , 
which means the potential surface describing the displacement of the lithium only could be 
either a single anharmonic well or a very shallow double well (corresponding to a positive or 
a negative sign). Linear and cubic terms (e.g; Q and Q 3 ) are excluded from the fit based on 
symmetry considerations, and terms which are not along the coordinates calculated; Qu, 
Qo and Qu,Oi are excluded from the fit since they degrade the variances of the quadratic 
coefficient. These include terms like Q\iQo an d QoQii- The dynamical matrix in units of 
Ryd 2 A~ 2 amu representing the LAPW energies of LiNbC>3 and LiTaC>3 are 



D LAPW (LiNb0 3 



' D Li D Li ,o^ 



( 



^ -0.001 -0.012^ 



/ 



V 



D LAPW (LiTa0 3 ) 



\ 



-0.012 -0.015 
\ 



0.0018 -0.012 
-0.012 -0.023 

The lithium-only contributions {Dn) are the smallest (an order of magnitude smaller then 
the rest), the oxygens only contributions are larger, with Dq in LiTaC>3 larger than LiNbCV 
Whether the origin of lithium and oxygen coupling is Coulombic can be determined by 
looking at the Madelung contribution to the dynamical matrices. The Madelung energies 
were calculated using experimental positions and lattice parameters and full ionic charges. 
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The second derivatives at zero displacements, which are the elements in the dynamical 
matrix, were calculated numerically and are shown in Table [TTI]. In the case of LiNbOs and 
LiTaOs the Madelung contributions to the dynamical matrix are 



D M A D (LiNb0 3 



Dm ad (LiTaO^ ) 



( 0.038 0.004 ^ 



i 0.004 -0.171 j 

( 0.04 0.003 ^ 
y 0.003 -0.162 j 



As expected, Do has the largest magnitude, followed by Du- The coupling term between 
the lithiums and the oxygens is in fact zero. This means that the origin of the lithium-oxygen 
coupling is not pure Coulombic (Madelung). 

Another possibility is the polarization of the oxygens by the lithium displacement, leading 
to changes in Nb(Ta)-0 bonding. The experimental ferroelectric configuration generates an 
effective dipole at the lithium sites. This dipole field can polarize the oxygens and drive 
them off center, yielding a ferroelectric distortion. We have plotted the self-consistent charge 
densities in two configurations; in one only the oxygens are displaced and in the other both 
the oxygens and the lithiums are displaced. In order to see the effects of displacing the 
lithiums we subtracted the two charge densities. This is shown in Fig. [5] where the charge 
density contours are plotted on a scale of -0.1 to 0.1 electrons/bohr 3 and the contour interval 
is 0.002 electrons/bohr 3 . A large dipole is seen at the lithium sites due to the displacement of 
the lithiums. Little polarization of the oxygens is observed; there is no evidence for any large 
dynamical covalency effects. We can therefore eliminate the possibility of oxygen-lithium 
coupling through either Madelung or polarization effects. 

Another possible source for the oxygen-lithium coupling is through the crystal structure. 
It is important to notice that the oxygens move not only along the c axis, but rather have 
sizable displacements along the a and b well. We have tested the importance of 

these displacements by moving only the oxygens along the c component of the experimental 
ferroelectric displacement (the polar axis) . The resulting energy curve was far shallower then 
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the energy surface which resulted from moving the oxygens only along the experimental soft 
mode coordinate. This is shown in Fig. ^ where the upper curve represents the displacements 
of the oxygens along the polar axis only and the lower, deeper well represents the total energy 
when displacing the oxygens along the experimental ferroelectric distortion (along a,b and 
c axes) . The reason for these big energy differences can be seen from Table |TJ which shows 
the Ta-O, Nb-0 and Li-0 bond lengths. The ionic radii of Li is about 0.6 A , that of Ta 
or Nb is about 0.6 A and the ionic radii of oxygen is about 1.4 A @, making the sum of 
each pair (Li-O, Nb-0 and Ta-O) about 2.0 A . When the oxygens are displaced only along 
the c axis, the oxygen-Nb (Ta) separation becomes only 1.83 (1.86) A which is about 0.17 
(0.14) A shorter then the sum of the ionic radii (about 8 (7)%). Therefore it is energetically 
favorable for the oxygens to displace in the a-b plane as they move along the c axis. 

If we now consider the experimental ferroelectric displacement of only the oxygens and 
the experimental ferroelectric coupled displacement of the oxygens and lithiums, displacing 
the oxygens only results in a Li-0 separation which is also shorter than the sum of their ionic 
radii (Table |T|). This explains why the wells associated with the oxygen displacements alone 
are shallower than those obtained with the displacement of both the lithiums and oxygens. 
The origin of the Li-0 coupling is therefore the fact that motion of the oxygens alone yields 
a Li-0 distance that is larger then the sum of their ionic radii, resulting in a deeper well for 
the coupled motion (in which the Li and oxygens move away from each other). 

We can therefore conclude that the driving mechanism behind the phase transformation 
in these systems is the displacement of the oxygens towards the transition metal atoms. 
Displacement of the oxygens in the direction of the transition metal atoms only (the c axis) 
would result in too short Nb (Ta)-oxygen bonds. The oxygens therefore move also in the 
plane perpendicular to the c axis, towards the lithiums. This shortens the lithium-oxygen 
bond so that the lithium displacements are coupled with the oxygen motions. 

The transition temperature, Tc, cannot be calculated directly from the zone center en- 
ergetics. In the usual models for ferroelectric phase transitions, Tq is related to the relative 
strength of the local (on-site) and coupling terms in the energy^. Since we find the zone 
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center energetics to be similar, the difference in Tc must be due to differences in the ener- 
getics at the zone boundary. 

In order to understand the origin of the ferroelectric distortion, we next examine the elec- 
tronic structure of these materials. 



B. Electronic Structure 

One goal of this research is to understand the origin of ferroelectricity in LiTaOs and 
LiNbOs and the difference in Tc from their electronic structure. Fig. |6]a compares the 
electronic density of states for LiTaOs and LiNbOs, both at the ferroelectric configuration. 
The energy scales are lined up with the top of the valence bands at zero energy. It is clear 
that the total density of states of these two materials is very similar. We can look further 
at the different contributions to the density of states; Fig. [|b compares the partial density 
of the Ta 5d state and the Nb 4d state both in the ferroelectric phase. The top of the 
valence band is composed mostly of oxygen p states. This figure shows a large density of 
transition metal d states in the valence band which means that the oxygen p states in these 
two materials are hybridized with the d states. The Nb d states have a large peak at the 
bottom of this band which is missing in the case of the 5d states of the Ta atom. The 
origin of this peak is the fact that the lowest valence bands of LiNbOs (the bands at about 
-4.5 eV or -0.35 Ryd) are less dispersive than the lowest valence bands of LiTaOs. This 
will be further discussed later. The same conclusion is derived from Fig. |]c of the partial 
density of the oxygen p state of the two materials in the valence band in the ferroelectric 
structure. Here too, the densities of states are very similar. The same peak at the bottom 
of the band is seen here for the LiNbOs, which is missing in the valence band of LiTaOs. 
All three figures that compare the total and partial density of states of the two materials in 
their ferroelectric phase show large hybridization between the transition metal d states and 
the oxygen p states, which is the reason for the oxygen displacements towards the transition 
metal atoms. 
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Next we compare the densities of states in the paraelectric and the ferroelectric phases. 
Fig. 0a illustrates the total density of states of LiNb0 3 in the paraelectric (solid line) and the 
ferroelectric (dashed line) phases. The bands in the two phases look similar except that the 
bands at the ferroelectric phase are slightly wider then the bands at the paraelectric phase. 
Fig. 0b compares the Nb 4d state both in the paraelectric (solid line) and the ferroelectric 
(dashed line) phases and Fig. 0c shows the Ta 5d states in the two configurations. The 
large peak at the lower part of this band (the peak at about -0.35 Ryd or -4.5 eV) is shifted 
in the ferroelectric case to higher energies. Fig. [7|d shows the density of oxygen p states of 
LiNbOs at the two phases. The peak at the bottom of the band is shifted in the ferroelectric 
phase from the paraelectric one. These figures indicate that the electronic structure at the 
paraelectric and the ferroelectric phases are quite similar. Fig.[7|e compares the lithium 2s 
character in the paraelectric and ferroelectric phases. It is evident that the lithium is almost 
completely ionized and that its electronic distribution does not change during the phase 
transition. 

Fig. |8| shows the band structure of LiNbOs in the ferroelectric state. The band gap is 
indirect, the top of the valence band is between T and Z, and the bottom of the conduction 
band is at the T point. The Brillouin zone for the rhombohedral lattice^ is illustrated in 
the inset in the figure. The energy between the Z and the A point was calculated along 
a straight line between the two points not along the Brillouin zone face for the purpose of 
comparison with Ching et a/.'s results!!. 

The band gap is 3.1 eV which is about 15 % lower than the value obtained from optical 
measurements of the near stoichiometric sample!!!! of 3.78 eV. The lithium 2s states are 
separated by 13.6 eV from the oxygen 2s states. These bands would not appear in the 
LiTaOa band structure since the lithium 2s states were treated as core states; these bands 
are very flat. The oxygen 2s are separated by 10 eV from the valence bands. The lowest 
conduction bands are the Nb 4d states for LiNbOs or Ta 5d states for LiTaOs. We have also 
compared this band structure with Ching et a/.!! who used the OLCAO method and got a 
band gap of 3.56 eV and the bands compare well with our results. 
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Fig. ^a and b show the band structure for LiTaOa and LiNb03, respectively. Each figure 
shows the ferroelectric (solid line) and the paraelectric (dashed line) phases. The changes 
observed between the ferroelectric and paraelectric bands are the band gap which is larger in 
the ferroelectric phase by about 15 % and the band width, especially the conduction band, 
which is larger in the paraelectric phase. 

The band structures of LiTaO"3 (dashed line) and LiNbO"3 (solid line) both in the ferro- 
electric phase, are shown in Fig. |H| The only difference between these two band structures 
is the larger band gap, by about 1 eV (30 %), in LiTaOs. The conduction bands are shifted 
by 1 eV from each other, but otherwise, their structure is almost the same. The valence 
bands are almost identical, which is consistent with the results of the total energy calcu- 
lations where the two well depths were found to be very close to each other (within 1.2 
mRy of each other) and the fact that the number of valence electrons in the muffin tins in 
both materials was similar. The only difference which was found between the two electronic 
structures was in the size of the band gap and the less dispersive nature of one band. The 
difference in the band gap will have an effect on quantities which include summing over 
unoccupied states as well as the occupied ones, like the polarizability. This difference will 
lead to different phonon dispersion in the two materials, and thus to different T^'s. Zone 
boundary or linear response calculations are necessary to further explore this issue. 



C. Comparison to the Perovskites 

Previously we studied the difference between the self-consistent charge densities and 
charge densities computed using overlapping ions with the PIB model for both LiTa03 
and LiNbOa in the ferroelectric phasell. In the PIB model which is a non-empirical ionic 
modefi, the charge densities are calculated via a Gordon-Kim type modelS, where the ions 
are allowed to breath corresponding to changes in the crystal potentials. The comparison 
indicated large hybridizations between the Ta atoms and the oxygens and between the Nb 
atoms and its oxygen neighbors and the Li atoms were fully ionized in the self-consistent 
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charge density. 

These results are consistent with the energetics and electronic structures results, all 
pointing to the same conclusion that the driving mechanism behind the ferroelectric insta- 
bility in the LiNb(Ta)03 systems is the hybridization between the d states on the transition 
metal atoms and the 2p states on the oxygens. The lithiums are but passive players in 
the ferroelectric instability. This is very similar to the ferroelectric mechanism in the per- 
ovskite ferroelectrics, where the oxygen-transition metal atom hybridization, in addition to 
the Coulombic long range interaction which tend to drive the system off center, overcome the 
short range repulsions which tend to leave the system in it's high-symmetry configuration. 

An interesting comparison can be made with the K(Nb,Ta)03 system; one major differ- 
ence is the fact that the perovskite KTa03 is an incipient ferroelectric where LiTaOa has a 
high transition temperature. In this sense a qualitative comparison can be made between the 
two sets of systems, as in both systems the transition temperature is higher in the niobate 
systems, being zero for KTa03. This would mean a shallower well for the tantalates, where 
in the case of KTa03 the well is apparently lower than the thermal vibrations, as shown by 
this study for the Li(Nb,Ta)03 systems and by Postnikov for the K(Nb,Ta)C>3 systems@. 

The electronic structures of the two sets of systems (LiNb(Ta)0 3 and KNb(Ta)0 3 ) show 
a large hybridization between the transition metal atoms and the oxygens, and the amount 
of hybridization between the transition metals and the oxygens in the two sets of systems 
is similar. This can be seen from Fig. which shows both the partial density of Nb 4d 
states and that of Ta 5d states in the valence bands in the ferroelectric configuration in 
the LiNb(Ta)0 3 systems. We can compare these results to the same densities of states 
calculated for KNb0 3 and KTa0 3 by LMTOfl reproduced in Fig. [TT| These figures show 
the Nb (Ta) density of states with the Nb (Ta) atom undisplaced and displaced by 0.073a (a 
being the lattice constant) along the <111> direction, which is an exaggerated displacement 
used to enhance the differences between the two phases. We can see that the same trends 
exists in this picture as in Fig. |6|b for the Nb (Ta) in LiNb03 (LiTaOs). The ferroelectric 
peaks at the bottom of this band are slightly shifted towards higher energies. In both cases 

16 



there is large hybridization between the transition metal atoms and the oxygens, but this 
hybridization does not change much during the transition. This hybridization is essential for 
the onset of the ferroelectric instability, however the amount of hybridization in this system 
doesn't change much through the phase transition like in other ferroelectrics, e.g BaTi03 or 
PbTi0 3 i. 

We find that LiNbOs and LiTaOs are almost identical in their electronic behavior. The 
amount of d character in the valence bands, which is a measure of the hybridization between 
the transition metal ions and their oxygen neighbors is very similar. This is in contrast to 
the conclusion that Ta is less ionic than Nb, reached by Postnikov et al. for K(Ta,Nb)C>3. 
This conclusion was based on the smaller transition metal atom density of states peak at 



the bottom of the valence band, observed in Figures |TT| a and b. It is seen from Fig. [T(] that 
the origin of this peak is the less dispersive nature of one LiTaC>3 band versus LiNbCV 

The similarity between LiNb(Ta)03 and KNb(Ta)3 is the fact that the driving mechanism 
for the phase transition in the two systems is oxygen-B atom hybridization. The difference 
between the two systems lies in the different structure which yields a different oxygen-A 
atom interaction. 

In both LiNb(Ta)03 and KNb(Ta)C>3 the hybridization between the B atoms (the tran- 
sition metals) and the oxygens causes the oxygens and the B atoms to displace towards 
each other. In the LiNb(Ta)03 system the oxygen-B atom separation is larger than the 
sum of their ionic radii and the oxygens markedly displace in the a-b plane as they move 
along the polar axis. This however, makes the oxygen-A atom separation larger then the 
sum of their ionic radii resulting in the coupled oxygen-A atom motion. This is in contrast 
to the perovskites where the A site is large enough to allow the oxygens to move towards 
the A atoms, e.g. KNbC>3 where the potassium-oxygen separation in the highest and lowest 
symmetry structures are about 2.85 A and 2.83A , respectively, compared with the sum of 
their ionic radii which is about 2.78 A . 

It is interesting to note that when doping KTa03 with lithium atoms (KLT), the system 
does displace off-center, with a critical concentration of lithiums as small as 2.2%B This 
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could be the result of the lithium ion having a much smaller ionic radii then the potassium 
with respect to the perovskite structure, being about 0.6 A for Li and 1.4A for kH. This 
would allow the lithiums, driven by Madelung forces, to displace off-center, and due to the 
large space open to the lithiums in the perovskite structure their amplitudes will be much 
larger then in the LiTaC>3 system, resulting in a dipole field that polarizes the oxygens and 
distorts them into off-center positions. The phase transition in KLT is significantly different 
then in a conventional ferroelectric and there is some discussion of whether KLT is a true 
ferroelectric. DiAntonio et a/.0 concluded that the coincidence of the temperature of the 
maximum of the dielectric permittivity with the appearance of other anomalies that are 
characteristics of a structural transformation are a sign that the transition is ferroelectric, 
whereas Azzini et al® state that the size of the domains having a homogeneous spontaneous 
polarization is significantly smaller than the size of the structural domains. 

In the Slater picture of the so called "rattling ion", the B atom lies off center because 
it is too small to fit into the oxygen octahedra surrounding it. This is in fact the opposite 
of the picture in the LiNb(Ta)C>3 systems where the separation of the oxygens from the Nb 
(Ta) atoms are smaller (1.9A) then the sum of their ionic radii (2A). Also, a comparison of 
the Ta-0 distances in KTaOs and LiTaC>3 shows exactly this same effect. In KTa03, the 
oxygen octahedra is larger than the oxygen octahedra in LiTaC>3, and yet, in KTa03 the 
B atom never displaces to the off center position, while in LiTaOs, the B cation exhibit a 
ferroelectric distortion. 

IV. CONCLUSIONS 

It is shown that LiNbC>3 and LiTaOs are very similar in both their electronic structure and 
energetics. The differences in their well depth are very small, the amount of hybridization in 
the two materials is similar and the charge densities are similar. Also, these two materials 
hardly change their electronic structure during a phase transition. The only difference found 
between these two systems is the difference in the conduction bands. Zone boundary effects 
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which are not included in this study and this difference in the electronic structure of the two 
systems are two possible candidates to explain the difference in the transition temperatures 
of the two systems. 

It is demonstrated that contrary to previous models which emphasized the hopping of 
the lithium atoms between the three positions as the driving mechanism for the phase 
transformation, in these systems, no triple well potential was found for the lithium motion. 
The deep double wells found are the result of the oxygen displacements towards the transition 
metal atoms, which are the result of the hybridization between the two atoms. The wells 
indicate an order-disorder character for the oxygen. Local changes in the oxygen octahedra 
are responsible for the lithium displacements from their centrosymmetric sites. The lithiums 
themselves are passive players in the ferroelectric energetics. 
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TABLES 

TABLE I. Total Energies for different configurations, a) LiNb03. b) LiTaC>3. Coordinates 
are in primitive rhombohedral coordinates. The first column refers to the displaced atom and the 
amount of displacement as a fraction of the paraelectric to experimental ferroelectric normal mode 
amplitude. 
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h.) LiTaO, 
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Li - 1.0 
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a Calculated using 28 k-points in the irreducible Brillouin zone to test convergence. Other points 
included 10 k-points. 

b This is the experimental ferroelectric disortion. 

c Energy shift due to different muffin tin sizes is included as described in the text. Energy shift is 
6.64 mRy. 

d Energy shift due to different muffin tin sizes is 6.6 mRy. 
e Energy shift due to different muffin tin sizes is 7.14 mRy. 
f Energy shift due to different muffin tin sizes is 9.68 mRy. 
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TABLE II. Bond length, in A, of transition metals-oxygens and Li-oxygens in different config- 
urations. Sum of ionic radii of each pair is about 2 A. 
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TABLE III. Parameter table for the polynomial fit of the LAPW and Madelung energies. 
Energies are in Ryd. LAPW LiNb0 3 energies are shifted by -16193 Ryd and LAPW LiTa0 3 
energies by -63395 Ryd. 



Coefficient 


LAPW energies 


Madelung energies 


Const 


-0.8467(3) 


-31.7268(0) 


Qh 


-0.0005(4) 


0.0194(0) 


Qh 


-0.0075(9) 


-0.0857(1) 


QuQo 


-0.0115(8) 


0.0040(0) 


QhQo 


0.0017(3) 




Qh 


0.0003(1) 




Qh 


0.0038(5) 






R 2 = 0.997 


R 2 = 0.999 


b.) LiTaCh 


Coefficient 


LAPW energies 


Madelung energies 


Const 


-0.0800(11) 


-31.7714(0) 


Qh 


0.0009(14) 


0.0200(8) 


Qh 


-0.0117(30) 


-0.0812(13) 


QhiQo 


-0.0121(30) 


0.0030(8) 


QhQo 


0.0010(11) 




Qh 


0.0001(3) 




Qh 


0.0090(2) 






R 2 = 0.968 


R 2 = 0.996 



27 



FIGURES 



FIG. 1. The a) paraelectric and b) ferroelectric structures of LiTa03 and L1NDO3. The hexag- 
onal unit cell is outlined. 



FIG. 2. a.) Potential energy surfaces of LiTaOs. The upper curves represent displacements of 
the Li atoms along the soft-mode coordinate, the middle shallow double wells represent displace- 
ments of the oxygens alone and the bottom curves represent the displacements of oxygens and Li 
atoms along the same coordinate. The curves represent a 4th order polynomial fit to the data. 
They were not constrained to go through the zero of energy. The abscissa represents displacement 
of the oxygen atoms from the paraelectric configuration, in A. b.) The same for LiNbOs. 



FIG. 3. The energy as a function of normal coordinate, fitted to a quadratic. The lines are the 
eigenstates for the 1-D independent harmonic oscillators. The difference between the ground state 
and the lowest excited state gives a frequency of 270 cmT x for LiTa03 and 250 cm -1 for LiNbC>3. 
Both are in very good agreement with experimental results, a) Energy versus normal coordinate 
for LiTa0 3 and b) The same for LiNb0 3 . 
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FIG. 4. Total energy surfaces of LiNbO-3 with only oxygens displaced along soft mode coordi- 
nate (lower curve) and with only the oxygens displaced along the c axis only (upper curve). The 
curves are a 4th order fit to the data. The abscissa represents displacements, in A, of the oxygens 
from their paraelectric positions. 



FIG. 5. L1NDO3: Charge density resulting from subtracting the charge density of a configura- 
tion in which oxygens only are displaced from the charge density of the full ferroelectric distortion 
(both lithiums and oxygens are displaced). The scale is from -0.1 to 0.1 electrons per bohr 3 and the 
contour interval is 0.002 electrons/bohr 3 . No evidence for dynamical covalency effects that would 
lead to coupling of oxygen and lithium motions are seen 



FIG. 6. a) Electronic density of LiTa03 (solid line) and LiNbOs (dashed line), both in the 
experimental ferroelectric configuration, b) Density of Ta 5d states (solid line) and Nb 4d states 
(dashed line) in the valence band, both the experimental ferroelectric configuration, c) Oxygen p 
states of LiTa03 (solid line) and of LiNbOs (dashed line) in the valence band, in the ferroelectric 
phase. 
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FIG. 7. a) Electronic density of states for LiNb03 in the paraelectric (solid line) and ferro- 
electric (dashed line) configurations, b) The Nb 4d state in the two phases, c) The LiTaOs Ta 
5d states in the two phases, d) The LiNbCVs p states on the oxygens in the paraelectric and the 
ferroelectric phase, e) The s orbital on the Li atom of LiNb03 in the two phases. The valence 
band top is lined up with the zero of energy. 



FIG. 8. The band structure of LiNb03 in the ferroelectric phase. The energy scale is in eV 
and the fermi level is shown. The band gap is 3.1 eV. The Li s states do not interact with the rest 
of the bands and are about 16 eV below the Oxygen 2s states. Inset: The brillouin zone. Some 
high symmetry points are illustrated. The energy between the Z and the A point was calculated 
along a straight line between the two points. From [60]. 



FIG. 9. a) The band structure of LiTa03 in the ferroelectric phase (solid line) and the para- 
electric (dashed line). Only the valence and the conduction bands are shown. The band gap in the 
ferroelectric phase is about 4.0 eV, and is decreased in the paraelectric phase by about 15 %. b) 
The same for LiNb03. The gap decreases by about 15 % between the ferroelectric and paraelectric 
phases. No major differences are observed in both systems between the two phases. 
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FIG. 10. The band structure of LiNbC>3 (solid line) and LiTaC>3 (dashed line), both in the 
ferroelectric structure. The LiTaC>3 band gap is larger than the LiNb03 band gap by about 1 eV. 
The valence bands of the two materials are almost identical. 



FIG. 11. From LMTO results of KTa0 3 and KNb0 3 , by Postnikov et al. [29]. a) Local 
densities of states at the Nb site with the Nb undisplaced (solid line) and displaced (dashed line) 
from it's rhombohedral position. The displacement of the Nb atom is exaggerated in order to 
enhance the trends shown, b) at the Ta site, under the same conditions. The units are in Ryd for 
the energies and Ryd -1 for the density of states. 
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